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Abstract. We speculate that a new generation of particle astrophysics instrumentation will 
reveal the enigmatic sources of cosmic rays prior to the one hundredth anniversary of their 
discovery by Hess in 1912. While only a "smoking gun" is missing for the case that the galactic 
component of the cosmic ray spectrum originates in supernova remnants, deciphering the origin 
of the extragalactic component is still at a level of reading tea leaves. 



1. The energetics of cosmic ray sources 

Cosmic accelerators produce particles with energies in excess of 10^ TeV; we still do not know 
where or how [1]. The flux of cosmic rays observed at Earth is shown in Fig. 1. The energy 
spectrum follows a broken power law. The two power laws are separated by a feature dubbed 
the "knee" at an energjj^ of approximately 3 PeV. There is evidence that cosmic rays up to 
this energy are Galactic in origin. Any association with our galaxy disappears in the vicinity 
of a second feature in the spectrum referred to as the "ankle"; see Fig. 1. Above the ankle, 
the gyroradius of a proton in the Galactic magnetic field exceeds the size of the Galaxy and 
it is routinely assumed that we are witnessing the onset of an extragalactic component in the 
spectrum that extends to energies beyond 100 EeV. Direct support for this assumption now 
comes from two experiments pi that have observed the telltale structure in the cosmic ray 
spectrum resulting from the absorption of the particle flux by the microwave background, the 
so-called Greissen-Zatsepin-Kuzmin cutofl^ The origin of the flux in the intermediate region 
covering PeV to EeV energies remains a complete mystery. Although the routine assumption 
is that it results from some high energy extension of the reach of the Galactic accelerators, no 
convincing mechanism for this has been identified. 

It is assumed that cosmic rays originate in cosmic accelerator^ Acceleration to TeV energy 
and above requires massive bulk fiows of relativistic charged particles. These are likely to 

^ Written in celebration of the 75th birthdays of Frank Avignone and Ettore Fiorini and in memory of the 75th 
birthday of the late Peter Rosen. 

^ We will use energy units TeV, PeV and EeV, increasing by factors of one thousand from GeV-energy. 

The possibility has not been eliminated that the "cutofF is simply the maximum energy reached by the 
accelerator(s). 

^ Speculations on the origin of the highest energy cosmic rays fall into two categories, top-down and bottom- 
up. In top-down models it is assumed that the cosmic rays are the decay products of cosmological remnants or 
topological defects associated with, for instance. Grand Unified theories with unification energy MQyp ~ 10^* eV. 
These models predict, besides cosmic rays, large fluxes of gamma rays and neutrinos that have not been observed. 
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Figure 1. At the energies of interest here, the cosmic ray spectrum follows a sequence of 3 power 
laws. The first two are separated by the knee, the second and third by the ankle. The cosmic 
rays beyond the ankle are a new population of particles produced in extragalactic sources. 



originate from the exceptional gravitational forces in the vicinity of black holes or neutron stars. 
Gravity powers large currents of charged particles that are the origin of high magnetic fields. 
These create the opportunity for particle acceleration by shocks. It is a fact that black holes 
accelerate electrons to high energy; astronomers observe them indirectly by their synchrotron 
radiation. Some must accelerate protons because we detect them as cosmic rays. 

The Swiss astronomer Fritz Zwicky suggested as early as 1933 that supernova remnants could 
be sources of the Galactic cosmic rays. It is assumed that the accelerators are powered by the 
conversion of 10^'' erg of energy into particle acceleration by diffusive shocks associated with 
young (~ 1000 year old) supernova remnants expanding into the interstellar medium f5]. Like a 
snowplough, the shock sweeps up the ~ 1 proton/cm^ density of hydrogen in the Galactic plane. 
The accumulation of dense filaments of particles in the outer reaches of the shock, clearly visible 
as sources of intense X-ray emission, are the sites of high magnetic fields. It is theorized that 
particles crossing these structures multiple times can be accelerated to high energies following 
an approximate power-law spectrum dN/dE ~ E^"^. The mechanism is familiar from solar flares 
where filaments of high magnetic fields accelerate nuclear particles to tens of GeV. The higher 
energies reached in cosmic ray accelerators are the consequence of particle flows of much larger 
intensity powered by the gravitational energy of collapsed objects such as neutron stars and 
black holes. 

From a myriad of ideas, speculations on the sites for the acceleration of extragalactic cosmic 



rays have converged on the supermassive black holes at the centers of active galactic nuclei 
(AGN) or the primary engines of gamma ray bursts (GRB). It should surprise nobody, however, 
if the final answer turned out to be neither of these. 

AGN are of special interest because some emit most of their luminosity at TeV energy and 
above. Their inferred isotropic luminosities can be as high as 10'^^ — lO^^ergs^^. They produce a 
typical spectrum dN/dE^ oc in the MeV-GeV range. The energetics requires mass accretion 
on a black hole that is up to one billion times more massive than our Sun. Some of the matter 
falling onto the black hole is deflected and accelerated in highly beamed jets aligned along 
its rotation axis. Both the inflow onto the central engine and the jet provide opportunities 
for the buildup of large magnetic fields and shock acceleration. For all these reasons AGN 
were pinpointed by Ginzberg and Syrovatskii j4j as far back as 1964 as candidate cosmic ray 
accelerators. A subset, called blazars, emit high-energy radiation in collimated jets pointing at 
the Earth and are the sources of photons with energies of tens of TeV. Their emission is highly 
variable over several time scales. TeV-energy bursts as short as minutes have been observed. 
Exceptionally, the nearest active galaxies. Centaur us A (Gen A) and M87 have been detected 
in TeV gamma rays even though their jets are not pointing at us. 

Some argue instead that gamma ray bursts (GRB), outshining the entire Universe for the 
duration of the burst, are the best motivated sources of high energy cosmic rays [5j. The collapse 
of a massive star to a black hole has emerged as the origin of "long" GRB with durations of 
tens of seconds. In the collapse a fireball is produced of electromagnetic plasma that expands 
with a highly relativistic velocity powered by radiation pressure. The fireball eventually runs 
into the stellar material that is still accreting onto the black hole. If it successfully punctures 
through this stellar envelope the fireball emerges to produce a GRB display. While the energy 
transferred to highly relativistic electrons is observed in the form of synchrotron radiation, it is 
a matter of speculation how much energy is transferred to protons. The assumption that GRB 
are the sources of the highest energy cosmic rays determines the energy of the baryons in the 
fireball. Accommodating the observed energy spectrum of extragalactic cosmic rays leads to the 
requirement of roughly equal efficiency for conversion of fireball energy into the kinetic energy 
of protons and electrons. 

It is routinely emphasized that the flux of cosmic rays, especially at the highest energies, is 
very low. For example, at the onset of the extragalactic component near lOEeV the flux is only 
at the level of one particle per kilometer squared per year for a typical array with a steradian 
acceptance in angle. This can be translated into an energy flux 

E {e'^\ = , ^"J""^ 7 ^ = 3 X 10-8 GeV cm-^ s-^ st'^ . (1) 

\ dE j (10i0cm2)(3 X 107s)sr ^' 

The particles may be few, but each carries an enormous energy that can be expressed in 
macroscopic units, tens of Joules for the highest energies. We can derive the average energy 
density Pe oi cosmic rays in the Universe using the relation that the total flux = velocity x 
density, or 
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We obtain 

47r /•^-- 3 X 10-\^ GeV iq TeV 

PE = — ^ dE ^-10-19^, (3) 

c JE^i^ E cur" cm-' 

taking the extreme energies of the accelerator(s) to be -E'maxZ-E'min — 10^. This is the 
corresponding value for cosmic rays of the energy density of, for instance, microwave photons of 
410 photons of 2.7 K per cubic centimeter]^ 

^ We note that the energy density of photons changes by less than 5 orders of magnitude over the electromagnetic 
spectrum, from radio waves to GeV-photons, while the flux drops by almost 20 orders of magnitude. 



The energy content derived "professionally" by integrating the observed extragalactic 
spectrum in Fig. 1, including the GZK feature, is ~ 3 x 10^^^ erg cm~^ [6 . This is within 
a factor of our back-of-the-envelope estimate (1 TeV = 1.6 erg). 

It has been realized for a long time that the corresponding quantity for the Galactic 
component of the spectrum may be a revealing number. If one repeats the integration that 
we just introduced for the Galactic flux in Fig. 1, from GeV-energy, below which it is shielded 
from us by the sun, up to the knee near 3PeV, one obtains an energy density of the cosmic rays 
in our Galaxy of 

Pe ~ 10~"'^^ergcm~^. 

This happens to be very close to the energy density of light in our Galaxy and to the energy 
density -B^/Svr in its microgauss magnetic field. This density, as well as the one for extragalactic 
cosmic rays in the Universe, will represent informative benchmarks for speculating on the sources. 

2. Multi-wavelength Astronomy: Cosmic rays, gamma rays and neutrinos 

Although the origin of cosmic rays remains a matter of speculation almost one century after 
their discovery, particle astrophysicists have developed ambitious instrumentation that probes 
their sources with increased sensitivity detecting cosmic rays, gamma rays and neutrinos. In 
the 3-prong attack on the cosmic ray problem, TeV-astronomy is by far the most mature [7]. 
A new generation of ground based air Cherenkov detectors has revealed plausible candidate 
cosmic ray sources that had been, surprisingly, not previously identified at other wavelengths. 
However, the basic hurdle to conclusively identify the observed TeV gamma rays as the decay 
products of pions produced by a cosmic ray accelerator has not been overcome. Synchrotron 
radiation by energetic electrons, routinely observed in non-thermal sources such as supernova 
remnants, AGN and GRB, cannot be excluded as their possible origin. Neutrinos from the 
charged pions accompanying pionic gamma rays can provide incontrovertible evidence for cosmic 
ray acceleration. The predicted fluxes are small and difficult to detect [8]. Let's quantify the 
problem. 

How many gamma rays and neutrinos are produced in association with the cosmic ray beam? 
For orientation, consider a neutrino beam produced at an accelerator laboratory; see Fig. 2. The 
accelerated protons interact with a target, referred to as a beam dump, producing pions that 
decay into gamma rays and neutrinos]^ Generically, a cosmic ray source should also be a beam 
dump. Cosmic rays accelerated in regions of high magnetic fields near black holes inevitably 
interact with the radiation surrounding it, for instance, UV photons in active galaxies or the 
MeV-photons in GRB fireballs; see Fig. 2. In these interactions they generate neutral and 
charged pions by the processes 

p + 7 ^ TT^ + p and p + 7 ^ — > vr^ + n. 

While the secondary protons may remain trapped in the high magnetic fields, neutrons and the 
decay products of neutral and charged pions escape. The energy escaping the source is therefore 
distributed among cosmic rays, gamma rays and neutrinos produced by the decay of neutrons, 
neutral pions and charged pions, respectively. 

In reality the relation between the secondary fluxes is more complex. The cosmic rays may 
interact with gas besides radiation producing equal numbers of pions of all three charges in 
hadronic collisions p + p n [ vr" + vr'^ + vr^] + X. Pionic gamma rays may cascade in the source 
steepening their spectrum and, as already pointed out, there may be an additional contribution 
to the gamma ray flux originating from purely electromagnetic processes. 

^ For earthbound accelerators the dump is designed to reabsorb all secondary electromagnetic and hadronic 
showers. Only neutrinos exit the dump. If Nature constructed such a "hidden source" in the heavens, conventional 
astronomy would not reveal it. 
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Figure 2. Sketch of a cosmic ray accelerator producing photons and neutrinos. Cosmic rays 
accelerated in the vicinity of a black hole produce pions in interactions with the radiation and 
gases surrounding it. Neutral and charged pions decay into photons and neutrinos, respectively. 

The secondary flux of neutrinos and gamma rays in astrophysical sources can be related to 
the initial accelerated proton spectrum by specifying the final-state multiplicities of the pions 
and their average fractional energies Xi relative to the proton. For instance, in hadronic collisions 
of cosmic rays interacting with gas in or near the source, equal numbers of each pion charge 
TT^ + vr"'" + 7r~ are produced. The neutral pions decay into two gamma-rays — > 7 -|- 7 and 
the charged pions into leptons and neutrinos vr^ — > -|- -|- -|- i^e. The effect of neutrino 
oscillations is to equalize the number of neutrinos of each flavor arriving at Earth, resulting into 
one neutrino of each flavor per charged pion. In summary, we obtain two-thirds of a charged 
pion per interacting proton and one neutrino (of each flavor) per charged pion; and one-third of 
a neutral pion per interacting proton and two photons per neutral pion. Therefore, 



The fractional energies Xi are determined by the average pion inelasticity measured to be 
approximately 0.2 at accelerators. In the approximation that the pion decay products carry 
equal energies, we obtain that x^, ~ 0.25 x — 0.05 and Xj ~ 0.5 x x^r ~ 0.1. 

Notice that one can eliminate the initial cosmic ray flux to obtain a relation between secondary 
gammas and neutrinos. It just follows from the fact that there are 2(1) charged pions for every 
neutral pion in pp(p7) interactions and, after taking into account oscillations, 2(8) photons for 
every i/^. The flux referred to sums over neutrinos and anti-neutrinos which the experiments 




(4) 



(5) 



cannot separate. For instance, for the case of a ii^ ^ spectrum elimination of the cosmic ray flux 
from above equations yields: 

In the photoproduction case the neutrino flux is further reduced by a factor 4. 
3. The Sources of Galactic Cosmic Rays 

We first recall that the energy density of the cosmic rays in our galaxy is pE ~ lO""*^^ ergcm"'^. 
Galactic cosmic rays are not forever; they diffuse within the microgauss fields and remain trapped 
for an average containment time of 3 x 10^ years. The power needed to maintain a steady energy 
density is therefore 10~^^ erg/cm'^s. For a nominal volume of the Galactic disk of 10^^ cm^ this 
requires an accelerator delivering lO^^erg/s. This happens to be 10% of the power produced 
by supernovae releasing 10^^ erg every 30 yearsj^ This coincidence is the basis for the idea that 
shocks produced by supernovae exploding into the interstellar medium are the accelerators of 
the Galactic cosmic rays. 

Despite valiant attempts, it has been impossible to conclusively pinpoint supernova remnants 
as the sources of cosmic rays by identifying accompanying gamma rays of pion origin. An 
alternative signature of comic ray acceleration may prove to be more powerful. Additionally, 
pions are produced when the cosmic rays in the expanding remnant interact with the hydrogen 
in the Galactic plane. In the star forming regions where supernova are more likely to occur, the 
cosmic rays can interact with dense molecular clouds. We next estimate the TeV radiation from 
the Galactic plane associated with this indirect signature tracing cosmic ray acceleration. 

A generic supernova remnant releasing an energy of ~ 10^*^ erg into the acceleration of 
cosmic rays will generate a source of TeV gamma rays by interacting with the hydrogen in the 
Galactic disk. The emissivity in pionic gamma rays is simply proportional to the density 
of cosmic rays Ucr and the target density n of hydrogen atoms. Here ricr — 4 x 10~^^cm~^ by 
integrating the spectrum in Fig. 1 for energies in excess of 1 TeV. For a E~'^ spectrum. 



Q7 - ^ { 1^ ) V"' ^cr (>1 TeV), (7) 



or 



~ 2cx^app [n Ucr] ■ (8) 
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The proportionality factor is determined by particle physics; x-y is the average energy of 
secondary photons relative to the cosmic ray protons and Xpp = (ncjpp)"^ is the proton interaction 
length {(Jpp ~ 40 mb) in a density n of hydrogen atoms. The corresponding luminosity is 

L^(>lTeV) ~ — ~ 10^^ photons s-\ (10) 

PE 

where W/ pe '^s the volume occupied by the supernova remnant. We here made the approximation 
that the volume of the young remnant is approximately given by W/pE or, that the density 
of particles in the remnant is not very different from the ambient energy density pE ~ 
lO^^^ergcm^^ of Galactic cosmic rays. 



^ 10^^ erg correspond to 1% of the binding energy of a neutron star after 99% is initially lost to neutrinos. 



We thus predict a rate of TeV photons from a supernova at a nominal distance d of order 
1 kpc of 



dN^ ^L^^ 12 _ 10-11 ( photons \ / W \ / / 

dE ^^^^'^''^ 47rd2 " V cm2s AlO^OergAlcm-sAlkpc; 

(11) 

Such sources must emerge in an all-sky TeV gamma ray survey performed with an instrument 
with the sensitivity of the Milagro experiment [9] . 

Above estimate ignores the fact that supernovae are not uniformly distributed but are 
associated with regions of star formation such as the Cygnus region and the center of the 
galaxy. Of particular interest is the relatively close association of O, B stars in the Cygnus 
region at only ~ Ikpc. While the density of matter in the galactic plane is one proton per 
cubic centimeter, dense molecular clouds are common in star forming regions and can reach 
densities that are larger by two or three orders of magnitude. These are therefore efficient at 
converting cosmic rays into pions that decay into gamma rays and neutrinos. They are an 
indirect signature for cosmic ray acceleration in nearby remnants and, unlike what is the case 
for the remnants themselves, there is no electromagnetic contribution to the TeV radiation that 
must be differentiated from the pionic gamma rays. The position of the knee in the spectrum 
indicates that some sources accelerate cosmic rays to energies of several PeV. These "PeVatrons" 
therefore produce pionic gamma rays whose spectrum extends to several hundred TeV without 
cutoff in interactions with the interstellar medium, in particular with the dense clouds |10j . 

To date, the Milagro collaboration has identified 6 such PeVatron candidates. Three of these 
sources, MGROJ1908+06, MGROJ2019+37, and MGROJ2031+41, have post-trial significances 
of > 4.5 a (the only other Milagro source of such statistical significance being the Crab nebula). 
Two of the remaining three hotspots, MGRO J2043-h36 and MGROJ2032-h37, are located within 
the Cygnus region of the Galaxy. The remaining hotspot, MGROJ1852-I-01, has a pre-trial 
significance of 5.1 cr; it is the strongest source but has the lowest significance, presumably because 
of large local background. One has to realize that, after subtracting the sources considered here, 
an excess of TeV gamma rays persists in the Milagro 's skymap from the general direction of the 
Cygnus region [11] . 

Of particular interest is MGROJ1908-I-06 which has been confirmed by the H.E.S.S. telescope. 
Their observations reveal a spectrum consistent with a E~'^ dependence from 500 GeV to 40 TeV 
without evidence for a cutoff [12]. In a follow-up analysis the MILAGRO collaboration |13| 
showed that its own data are consistent with an extension of the H.E.S.S. spectrum to at least 
90 TeV; see Fig. 3. This is suggestive of pionic gamma rays from a PeVatron whose cosmic 
ray beam extends to the knee in the cosmic ray spectrum at PeV energies. Another source, 
MGROJ2031-I-4, has been observed |14j by the MAGIC telescope with a spectrum that is also 
consistent with E~'^. The lower fiux measured by MAGIC can be attributed to the problem of 
differentiating the source from the background in a high density environment like the Cygnus 
region. 

In the end, despite the suggestive evidence, conclusively tracing the observed gamma rays to 
pions produced by cosmic-ray accelerators has so far been elusive. It is one of the main missions of 
neutrino telescopes to produce the smoking gun for cosmic-ray production by detecting neutrinos 
associated with the charged pions. Neutrino telescopes detect the Cherenkov radiation from 
secondary particles produced by the interactions of high energy neutrinos in highly transparent 
and well shielded deep water or ice. They take advantage of the large cross section of high-energy 
neutrinos and the long range of the muons produced. The instrumentation consists of a lattice 
of photomultipliers deployed in a shielded and optically clear medium that is transformed into a 
Cherenkov detector. The IceCube telescope fTS] is under construction and is taking data with 
a partial array of 2400 ten inch photomultipliers positioned between 1500 and 2500 meter and 




Figure 3. The gamma ray and neutrino fluxes from MGROJ1908+06. The shading surrounding 
the neutrino flux represents the range in spectra consistent with the statistical and systematic 
uncertainties on the input gamma ray flux. Also shown is the flux of atmospheric neutrinos 
at the same zenith angle as the source (dashed line), taking into account the source size and 
angular resolution. 



deployed as beads on 40 strings below the geographic South Pole. With the completion of the 
detector by 2010-2011 the volume presently instrumented will be doubled from 0.5 to Ikm^. 

As previously emphasized, particle physics is sufficient to compute the neutrino fluxes 
associated with pionic gamma rays. The neutrino flux associated with MGRO J1908-I-06 is 
shown in Fig. 3 along with the TeV gamma ray observations from which it is derived. The 
sensitivity of IceCube to the Milagro sources has been evaluated assuming that the 6 sources 
represent the imprint of the Galactic cosmic-ray accelerators on the TeV sky [16]. While the 
number of events with energies of tens of TeV is relatively low, it is clear that this is the 
energy region where the atmospheric neutrino background, also shown in Fig. 3, is suppressed 
and an excess from these sources can be statistically established. While observing individual 
sources may in some cases be challenging, establishing a correlation between the Milagro and 
IceCube sky maps should be conclusive after several years; see Fig. 4. A "stacked" source 
search that will look for correlations between all six Milagro sources and the IceCube sky map 
shown has a Poisson probability of 3-5 a after 5 years; for a sample calculation see Fig. 5. The 
range reflects the imprecise knowledge of the gamma ray fluxes. The use of optimised methods 
using unbinned searches beyond the simple binned method considered here will further increase 
IceCubes sensitivity. 

It is important to emphasize that the photon flux from the Milagro sources is consistent 
with the flux expected from a typical cosmic-ray-generating supernova remnant interacting with 
the interstellar medium. In other words, the TeV flux is consistent with the energetics that is 
required to power the cosmic ray flux in the Galaxy. Alternative sources such as microquasars 
have been theorized to contribute to the Galactic cosmic rays. If this were indeed the case, 
cosmic ray energetics would require that they leave their imprint on the Milagro skymap, but 
none have been observed so far. 
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Figure 4. Simulated sky map of IceCube in Galactic coordinates after 5 years of operation 
of the completed detector. Two of the Milagro sources are visible "by eye" with 4 events for 
MGROJ1852+01 and 3 events for MGROJ1908+06 with energy in excess of 40 TeV. These as 
well as the background events have been randomly distributed according to the resolution of the 
detector and the size of the sources. 




Figure 5. Poisson probability of the excess from the 6 Milagro hotspots after 5 years as a 
function of the energy of the neutrino-induced muons and a gamma-ray cut-off of the source at 
SOOTeV. Clearly most sensitivity is associated with muons with energies in excess of 40 TeV. 



4. The Sources of the Extragalactic Cosmic Rays 

We first return to the observable that seemed so revealing in the case of the Galactic component, 
the energy density which is ~ 3 x 10~^^erg cm~^ for extragalactic cosmic rays. The power 
required for a population of sources to generate this energy density over the Hubble time of 
10^*^ years is ~ 3 x 10^'' erg s~^ per (Mpc)^ or, as often quoted in the literature, ~ 5 x 10^^ TeV 
per (Mpc)^ per year. 

As previously discussed, a GRB fireball converts a fraction of a solar mass into the acceleration 
of electrons, seen as synchrotron photons. The energy in extragalactic cosmic rays can be 
accommodated with the reasonable assumption that shocks in the expanding GRB fireball 
convert roughly equal energies into the acceleration of electrons and cosmic rays. It so happens 
that ~ 2 X 10^^ erg per cosmological gamma ray burst will yield the observed energy density in 
cosmic rays after 10^*^ years given that their rate is of order 300 per Gpc^ per year. Hundreds of 
bursts per year over Hubble time produce the observed cosmic ray density, just like 3 supernova 
per century accommodate the steady flux in the Galaxy. Problem solved? Not really, it turns 
out that the same result can be achieved with active galaxies. 

The energy density of 3 x lO^-'^^erg cm~'^ works out to not only [6j 

• ~ 2 X 10^^ erg per cosmological gamma ray burst, but also to 

• ~ 3 X 10^^ erg s~^ per cluster of galaxies, 

• 2 X 10^^ erg s^^ per active galaxy. 

The coincidence between above numbers and the observed output in electromagnetic energy 
of these sources explains why they have emerged as the leading candidates for the cosmic ray 
accelerators. Whether GRB or AGN, the observational evidence that the sources radiate similar 
energies in photons and cosmic rays is consistent with the beam dump scenario previously 
discussed. In the interaction of cosmic rays with radiation and gases near the black hole, roughly 
equal energy goes into the secondary neutrons, neutral and charged pions whose energy ends up 
in cosmic rays, gamma rays and neutrinos, respectively. 

Can IceCube reveal the extragalactic cosmic ray sources? Naively, the neutrino flux should be 
the same as the observed flux of cosmic rays of Eq. 1 . Not so naively, it is about 5 times smallei|^ 
we will refer to the energetics estimate for the neutrino flux accompanying the extragalactic 
cosmic ray accelerators as the band shown in Fig. 6, or 



As illustrated in the figure, after 7 years of operation AMANDA's sensitivity is approaching the 
interesting range, but it takes IceCube to explore it. 

Recently the Pierre Auger observatory has presented evidence that the arrival directions of 
cosmic rays with energies in excess of 6 X 10^ TeV may be correlated with nearby active galactic 



nuclei p!9]. The skymap suggests that the events cluster in the direction of the supergalactic 



plana^ possibly in the direction of the FRI radio galaxy Centaurus A, by far the closest AGN 



at only 3.4 Mpc. If AGN are indeed the sources, Cen A's proximity singles it out as a primary 
accelerator. We will review the information on its spectral energy distibution and argue that the 
observations are consistent, within large ambiguities, with the neutrino flux estimated above. 
In fact, Cen A and M87, the next nearest FRI whose jet is not aligned to our line of sight, have 
previously been singled out as potential cosmic ray accelerators on the basis of gamma-ray data 



* For experts, it is the Waxman-Bahcall "bound" adjusted downward because only 20% of the proton energy is 
transferred to pions [17'. 

^ One should of course realize that any sources limited to a GZK absorption length of less than 100 Mpc are 
likely to be clustered in the supergalactic plane. 



EldN/dEy = 1 - 5 X 10"*^ GeV cm^^ s~^ sr 



-1 



(12) 





IceCube 



AMANDA 



evolution of sources 
- FRI 



3 



6 



9 

log(^/GeV) 



12 



15 



Figure 6. Our energetics estimate of the flux of neutrinos associated with the sources of the 
highest energy cosmic rays (the shaded range) is compared to the limits established by the 
AMANDA experiment and the sensitivity of IceCube fT8| . Also shown is the flux derived from 
the assumption that AGN are the sources which we model using spectral energy information on 
the nearby active galacties Cen A and M87; see Fig. 7. Integration of AGN to larger redshifts 
can reconcile the two estimates which differ by a factor 3. Also shown is the background flux of 
atmospheric neutrinos. 

The problem of interpreting TeV observations is however challenging because the high energy 
emission of AGN is highly variable and it is difficult to compare multi-wavelength data taken at 
different times. Extragalactic TeV sources observed so far are BLLacs, a subclass of Fanaroff- 
Riley I (FRI) active galaxies whose jets are oriented along our line of sight. The data have yielded 
no hints of cosmic ray acceleration so far; the observed spectra can be readily accommodated 
by synchrotron radiation by electrons, followed by inverse Compton scattering of the photons 
to TeV energy. Exceptionally, the nearby (16 Mpc) FRI M87 was observed in the late 1990's by 
the HEGRA stereoscopic system of five imaging atmospheric Cherenkov telescopes despite the 
fact that its jet is angled to more than 30° from our line of sight [21 . 

Where Cen A is concerned, there exists archival data of TeV emission collected in the early 
1970's with the Narrabri optical intensity interferometer of the University of Sydney p2j. At the 
time the Sydney interferometer unsuccessfully searched the sky for gamma ray sources detecting 
the highest fluctuation in the direction of Cen A. As a followup they exclusively observed the 
source for a 3 year period. A variable flux was observed in two periods of heightened activity 
lasting roughly one year, pointing at a region of coherent emission of size of order 0.3 pc. This is 
consistent with the idea that the high energy emission is from an isotropic region near the base 
of the jet or the central black hole of mass 2x10^ solar masses, about two orders of magnitude 
more massive than the one at the center of our galaxy |20j . The fact that the interferometer 
beam did not include the radio lobes at the end of the jets further supports the idea of a central 
engine at the base of the jets. An average flux of 



dN^/d{lnE){E^ > 0.3 TeV) = 4.4 ± 1.0 x lO^^cm'^s" 
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(13) 
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Figure 7. Spectral energy distribution of Cen A (black dots). Having in mind that the source is 
variable, we show our estimates for the flux of TeV gamma rays (gray shading) and cosmic rays 
assuming that between 1 and 10 events observed by Auger originated at Cen A (blue shading). 
We note that the cosmic ray and TeV gamma ray fluxes estimated in this paper are at the level 
of the electromagnetic component shown from radio waves to GeV photons. Our estimate for 
the neutrino flux (labeled upper limit) is shown as the red line. 

was reported for the period of 3 years over which the source was variable. 

The integrated flux of Cen A is close to the one observed by HEGRA [23j from M87 

dN^{E^/d(lnE) > 0.73 TeV) = 0.96 ± 0.23 x 10"^2cm"^s"\ (14) 

after scaling the flux of M87 at 16Mpc to the distance to Cen A and adjusting for the different 
thresholds of the experiments. The observations are therefore consistent with identical source 
luminosities of roughly 7 x lO^'^ergs"^, assuming an E~'^ gamma-ray spectrum. This suggests 
that they may be generic FRI, a fact we will exploit to construct the diffuse neutrino flux from 
all FRI. Recently, MAGIC reported short day-long bursts of M87 [21] pointing at acceleration of 
particles even closer to the black hole. This, as well as the absence of associated X-ray activity 
expected in the case of electromagnetic processes, further supports the possibility of cosmic ray 
acceleration. 

For completeness, the H.E.S.S. experiment obtained a limit on Cen A |25| of 10~^^ in the 
same units, a flux smaller than the one just argued for. Given the burst nature of the data the 
disagreement may be acceptable. We have compiled the TeV gamma ray data in Fig. 7 along 
with other observations on the multi-wavelength emission spectrum of Cen A recently compiled 
by Lipari p6|. 

The same conversion of TeV gamma rays to neutrinos, exploited for Galactic sources, yields 



a neutrino flux of 

^ < 5 X 10-13 (^^y\eV-'cu.-'s-\ (15) 

conservatively normalized not to exceed the contemporary limit of the H.E.S.S. experiment. 

Finally, having previously argued for a relation between the energetics of TeV gamma rays, 
cosmic rays and neutrinos, we ask the question what Auger data may reveal about the flux of 
the source. The answer depends on the number of events in the sky map, A'events) that actually 
correlate with Cen A, a number that depends on the angular broadening of the source by the 
deflection of the particles by magnetic fields. A'^events is therefore a matter of speculation. We 
estimate the flux assuming a power-law spectrum of the form dN/dE = Nq{E / Eq)°^ , where 
the normalization Nq is flxed by the number of events observed Aovents- For a = —2.0, 
we have Aevents = fleld of view x time x efficiency x AoZ-Ethrcsh- This gives us dNcr/dE = 
Aevents X lO^"'^'^ (ii^/TeV)~^ TeV-"*^ cm^^ s"^ . The number of events with energy in excess of 
6 X 10'' TeV may be as many as 10, thus obtaining the range of cosmic ray fluxes shown in Fig. 7. 
It is close to, but below the flux of TeV gamma rays and suggest an even smaller neutrino 
flux than the one derived from TeV gamma ray information. The variability of the source and 
the possibility of a more complicated shape of the spectrum clearly prevent us from reaching 
quantitative results on the basis of the present data. 

The neutrino flux from a single source such as Cen A is clearly small: repeating the calculation 
for power-law spectra between 2.0 and 3.0, we obtain, in a generic neutrino detector of effective 
muon area Ikm^, between 0.8 and 0.02 events/year only. Having estimated the neutrino flux 
from a point cosmic ray source, we calculate next the total diffuse flux from all such sources 
within our horizon. Given an FRI density of n ~ 8 x 10^ Gpc-^ within a horizon of i? ~ 3 Gpc 
|27j . the total diffuse flux from all 47rsr of the sky is simply the sum of the luminosities of the 
sources weighted by their distance: 

^ = yJ^ = L,nR = And^R^, (16) 

where is given by Eq. 1 . We performed the sum by assuming that the galaxies are uniformly 
distributed. This evaluates to: 

'^^'^ -2 X 10-9 ( T^V' GeV-^cm-^s-isr-i, (17) 



dEdis VGeV 

approximately a factor 3 below the flux estimated on the basis of source energetics. Varying 
the spectral indices as before, we obtain an event rate per km^ year from the northern sky of 
between 19 and 0.5 neutrinos per year. Considering sources out to 3 Gpc, or a redshift of order 
0.5 only, is probably conservative. Extending the sources beyond z ~ 1 increases the flux by a 
factor 3 or so as was discussed in connection with the results in Fig. 6. We can thus bridge the 
gap between this and the previous estimate of the cosmic neutrino flux based on the energetics 
of the sources. 

The predicted flux of Eq. 16 should be within reach of IceCube and a future Mediterranean 
kilometer-scale neutrino telescope. The flux is close to IceCube's 90% confidence level sensitivity 
in a single year. It follows that detection at the 5a level should occur within ~ 5 years 

We here used the atmospheric background flux from charm particles predicted by ,28; rather than the larger flux 
conservatively implemented in the sensitivity calculation presented in the IceCube publication. The enhancement 
in the charm flux relative to the one used here is a direct consequence of the assumption that charm fragmentation 
functions scale and therefore clearly unphysical. 



It is not challenging to produce models that generate a neutrino flux at the level proposed 
here. For instance [2D], the central engine in FRI such as Cen A and M87 may feed a beam 
dump consisting of the gas surrounding the supermassive black hol(|5 

In summary, while the road to the identification of the sources of the Galactic cosmic ray has 
been mapped, the origin of the extragalactic component remains as enigmatic as ever. 
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